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■ linearly independent pure quantum states. We show that the design of the optimal measurement 

that minimizes the probability of an inconclusive result can be formulated as a semidefinite 
CN| ' programming problem. Based on this formulation, we develop a set of necessary and sufficient 

conditions for an optimal quantum measurement. We show that the optimal measurement can 
be computed very efficiently in polynomial time by exploiting the many well-known algorithms 
Oh! for solving semidefinite programs, which are guaranteed to converge to the global optimum. 

Using the general conditions for optimality, we derive necessary and sufficient conditions 
so that the measurement that results in an equal probability of an inconclusive result for each 
one of the quantum states is optimal. We refer to this measurement as the equal-probability 
q-( measurement (EPM). We then show that for any state set, the prior probabilities of the states 

can be chosen such that the EPM is optimal. 

Finally, we consider state sets with strong symmetry properties and equal prior probabilities 
for which the EPM is optimal. We first consider geometrically uniform state sets that are 
defined over a group of unitary matrices and are generated by a single generating vector. We 
then consider compound geometrically uniform state sets which are generated by a group of 
unitary matrices using multiple generating vectors, where the generating vectors satisfy a certain 
(weighted) norm constraint. 



In this paper we consider the problem of unambiguous discrimination between a set of 
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1 Introduction 



In recent years, research into the foundations of quantum physics has led to the emerging field of 
quantum information theory Q . Quantum information theory refers to the distinctive information 
processing properties of quantum systems, which arise when information is stored in or retrieved 
from quantum states. To convey information using quantum states, we may prepare a quantum 
system in a pure quantum state, drawn from a collection of known states < i < m}. To 

detect the information, the system is subjected to a quantum measurement. If the given states 
\<pi) are not orthogonal, then no measurement can distinguish perfectly between them [||]. A 
fundamental problem therefore is to design measurements optimized to distinguish between pure 
nonorthogonal quantum states. 

We may formulate this problem within the framework of quantum detection, and seek the 
measurement that minimizes the probability of a detection error, or more generally, the Bayes cost 
||, ||, [|, b]. More recently, a different approach to the problem has emerged, which in some cases 
may be more useful. This approach, referred to as unambiguous discrimination of quantum states, 
combines error free discrimination with a certain fraction of inconclusive results. The basic idea, 
pioneered by Ivanovic [jffl, is to design a measurement that with a certain probability returns an 
inconclusive result, but such that if the measurement returns an answer, then the answer is correct 
with probability 1. Given an ensemble consisting of m states, the measurement therefore consists 
of m + 1 measurement operators corresponding tom + 1 outcomes, where m outcomes correspond 
to detection of each of the states and the additional outcome corresponds to an inconclusive result. 

Ivanovic j7j developed a measurement which discriminates unambiguously between a pair of 
nonorthogonal pure states. The measurement gives the smallest possible probability of obtaining 
an inconclusive result for unambiguous discrimination, when distinguishing between two linearly 
independent nonorthogonal states with equal prior probabilities. This measurement was then fur- 
ther investigated by Dieks S and Peres M, and was later extended by Jaeger and Shimony |1C] to 



the case in which the two states have unequal prior probabilities. 

Although the two-state problem is well developed, the problem of unambiguous discrimination 
between multiple quantum states has received considerably less attention. In Peres and Terno 



consider unambiguous discrimination between 3 quantum states. Chefles [12] showed that a nec- 
essary and sufficient condition for the existence of unambiguous measurements for distinguishing 
between m quantum states is that the states are linearly independent. He also proposed a simple 
suboptimal measurement for unambiguous discrimination for which the probability of an incon- 
clusive result is the same regardless of the state of the system. Equivalently, the measurement 
yields an equal probability of correctly detecting each one of the ensemble states. We refer to such 
a measurement as an equal-probability measurement (EPM). Chefles and Barnett fi"3fl developed 
the optimal measurement for the special case in which the state vectors form a cyclic set, i.e., the 
vectors are generated by a cyclic group of unitary matrices using a single generating vector, and 
showed that it coincides with the EPM. In their paper, they raise the question of whether or not 
this is the only case for which the EPM is optimal. 

In this paper we develop a general framework for unambiguous state discrimination which 
can be applied to any number of states with arbitrary prior probabilities. For our measurement 
we consider general positive operator- valued measures fi O], consisting of m + 1 measurement 
operators. We derive a set of necessary and sufficient conditions for an optimal measurement 
that minimizes the probability of an inconclusive result, by exploiting principles of duality theory 
in vector space optimization. In analogy to the quantum detection problem, deriving a closed- 
form analytical expression for the optimal measurement directly from these conditions is a difficult 
problem. However, our formulation has several advantages. First, it readily lends itself to efficient 
computational methods. Specifically, we show that the optimal measurement can be found by 
solving a standard semidefinite program (SDP) |15[| , which is a convex optimization problem. By 
exploiting the many well-known algorithms for solving SDPs [jl6|, |l7| , the optimal measurement can 
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be computed very efficiently in polynomial time. Since an SDP is convex, it does not suffer from 
local optimums, so that SDP-based algorithms are guaranteed to converge to the global optimum. 
Second, although the necessary and sufficient conditions are hard to solve directly, they can be 
used to verify a solution. Finally, the necessary and sufficient conditions lead to further insight into 
the optimal measurement. In particular, using these conditions we derive necessary and sufficient 
conditions on the state vectors, so that the EPM minimizes the probability of an inconclusive result. 
In contrast with the general optimality conditions, these conditions can be easily verified given the 
state ensemble and the prior probabilities. Using these conditions we show that for any set of state 
vectors the prior probabilities can be chosen such that the EPM is optimal. 

Based on the necessary and sufficient conditions we develop the optimal measurement for state 
sets with broad symmetry properties. In particular, we consider geometrically uniform (GU) state 



sets [18, 19, ^] defined over a group of unitary matrices. For such state sets we show that the 
optimal measurement is the EPM, and we obtain a convenient characterization of the EPM that 



exploits the state symmetries. We then consider compound GU (CGU) state sets [21, 20 1 in which 
the state vectors are generated by a group of unitary matrices using multiple generating vectors. We 
obtain a convenient characterization of the EPM in this case, and show that when the generating 
vectors satisfy a certain constraint, the EPM is optimal. 

The paper is organized as follows. After a statement of the problem in Section ||, in Section |3] 
we derive the necessary and sufficient conditions for the optimal measurement that minimizes the 
probability of an inconclusive result, by formulating the problem as an SDP. In Section |] we consider 
the EPM and derive necessary and sufficient conditions on the state set and the prior probabilities 
so that the EPM is optimal. Efficient iterative algorithms for minimizing the probability of an 
inconclusive result which are guaranteed to converge to the global optimum are considered in 
Section |5[ In Sections || and ^ we derive the optimal measurement for state sets with certain 
symmetry properties, and show that the optimal measurement coincides with the EPM. 
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2 Unambiguous Discrimination of Quantum States 

Assume that a quantum system is prepared in a pure quantum state drawn from a collection of 
given states < i < m} in an r-dimensional complex Hilbert space 7i, with r > m. The 

states span a subspace U of 7i. To detect the state of the system a measurement is constructed 
comprising m + 1 measurement operators {IT,0 <i< m} that satisfy 

m 

I> = I r . (1) 

i=0 

The measurement operators are constructed so that either the state is correctly detected, or the 
measurement returns an inconclusive result. Thus, each of the operators IT, 1 < i < m correspond 
to detection of the corresponding states \4>i),l < i < m, and Ho corresponds to an inconclusive 
result. 

Given that the state of the system is \4>i), the probability of obtaining outcome k is 
Therefore, to ensure that each state is either correctly detected or an inconclusive result is obtained, 
we must have 

((f>i\Ilk\<t>i) =pAki l<i,k<m, (2) 

for some < Pi < 1. Since from ([!]), Ho = I r — Y^ILi ^t? @ implies that (0j|no|<^i) = 1 — Pi, so 
that given that the state of the system is \4>i), the state is correctly detected with probability pi, 
and an inconclusive result is returned with probability 1 — pi . 



It was shown in [12| that (g) can be satisfied if and only if the vectors \<pi) are linearly indepen- 
dent, or equivalently, dimU = m. We therefore make this assumption throughout the paper. In 
this case, we may choose 

IT = Pi\4>i)(4>i\ = PiQi, l<i<m, (3) 
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where 

Qi = \4>i)(4>i\, l<i<m, (4) 

and the vectors \<pi) G hi are the reciprocal states associated with the states \4>i) i.e., there are the 
unique vectors in U such that 

{<t>i\<t>k) = $ik, l<i,k<m. (5) 
With $ and <3? denoting the matrices of columns \4>i) and \4>i) respectively, 

$ = $($*$)- 1 . (6) 

Since the vectors \(f>i) are linearly independent, $*$ is always invertible. Alternatively, 

i = ($$*)t$, (7) 

so that 

|&> = (<M>*)t|&>, (8) 



where {■)' denotes the Moore-Penrose pseudo-inverse [22|; the inverse is taken on the subspace 
spanned by the columns of the matrix. 

We can immediately verify that the measurement operators given by (||) satisfy If r = m 
so that the dimension of Ti is equal to the dimension of the space U spanned by the vectors \4>i), 
then these operators are the unique operators satisfying @. If on the other hand r > m, then the 
measurement operators are not strictly unique. Indeed, any measurement operators of the form 

^k=PiQi + \fJ>i)(m\i l<i<m, (9) 
where G hi , also satisfy (0). Since \<fii) G Z^, (4>i\^k) = for every i, fc so that the measurement 
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operators given by (||) and (|9|) lead to the same detection probabilities (4>i\Ilk\4>i) = Pidik- We may 
therefore assume without loss of generality that the operators II j are restricted to U, so that they 
have the form given by (H|). 

If the state \4>i) is prepared with prior probability r]i, then the total probability of correctly 
detecting the state is 

m m 

Pd = Y, ^(^|n#i) = Wi- (io) 

i=l i=l 

Our problem therefore is to choose the measurement operators Ilj = piQ u or equivalently the 
probabilities pi > 0, to maximize Pjj, subject to the constraint (||). We can express this constraint 
directly in terms of the probabilities pi as 

m m 

j2u i = j^ Pi Q i <i r . ( n ) 

i=l i=l 

Note that (|ll|) implies that pi < 1. 

3 Semidefinite Programming Formulation 

We now show that our maximization problem (|l0| ) and (|ll]) can be formulated as a standard 
semidefinite program (SDP) |l5|, O], which is a convex optimization problem. There are several 
advantages to this formulation. First, the SDP formulation readily lends itself to efficient computa- 
tional methods. Specifically, by exploiting the many well known algorithms for solving SDPs [O], 
e.g., interior point methods^ [pi], 0], the optimal measurement can be computed very efficiently 
in polynomial time. Furthermore, SDP-based algorithms are guaranteed to converge to the global 
optimum. Second, by exploiting principles of duality theory in vector space optimization, the SDP 

interior point methods are iterative algorithms that terminate once a pre-specified accuracy has been reached. 
A worst-case analysis of interior point methods shows that the effort required to solve an SDP to a given accuracy 
grows no faster than a polynomial of the problem size. In practice, the algorithms behave much better than predicted 
by the worst case analysis, and in fact in many cases the number of iterations is almost constant in the size of the 
problem. 
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formulation can be used to derive a set of necessary and sufficient conditions for the probabilities 
Pi to maximize Po of (|To| ) subject to the constraint (p"l|). 

We note that recently SDP based methods have been employed in a variety of different problems 



in quantum detection and quantum information M, 23, 24, BR 26, 27] 



After a description of the general SDP problem in Section , in Section ^2 we show that our 
maximization problem can be formulated as an SDP. Based on this formulation, we derive a set of 
necessary and sufficient conditions on the measurement operators, or equivalently, the probabilities 
Pi, to minimize the probability of an inconclusive result. Although in general obtaining a closed form 
analytical solution directly from these conditions is a difficult problem, the conditions can be used 
to verify whether or not a set of measurement operators is optimal. Furthermore, these conditions 
lead to further insight into the optimal measurement operators. In particular, in Section || we use 
these conditions to develop necessary and sufficient conditions on the state vectors and the prior 
probabilities so that the EPM is optimal. 

3.1 Semidefinite Programming 

A standard SDP is the problem of minimizing 

P(x) = (c\x) (12) 

subject to 

F{x) > 0, (13) 

where 

m 

F{x) = F + Y J *iFi- (14) 

i=l 
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Here \x) S TZ m is the vector to be optimized, Xi denotes the «th component of \x), |c) is a given 
vector in lZ m , and F{ are given matrices in the space B n of n x n Hermitian matrices^]. 

The problem of (|l^) and ( |l3|) is referred to as the primal problem. A vector \x) is said to be 
primal feasible if F{x) > 0, and is strictly primal feasible if F(x) > 0. If there exists a strictly 
feasible point, then the primal problem is said to be strictly feasible. We denote the optimal value 
of P(x) by P. 

An SDP is a convex optimization problem and can be solved very efficiently. Furthermore, 
iterative algorithms for solving SDPs are guaranteed to converge to the global minimum. The 
SDP formulation can also be used to derive necessary and sufficient conditions for optimality by 
exploiting principles of duality theory. The essential idea is to formulate a dual problem of the form 
maxz D(Z) for some linear functional D whose maximal value D serves as a certificate for P. That 
is, for all feasible values of Z £ B n , i.e., values of Z £ B n that satisfy a certain set of constraints, 
and for all feasible values of \x), D(Z) < P(x), so that the dual problem provides a lower bound 
on the optimal value of the original (primal) problem. If in addition we can establish that P = D, 
then this equality can be used to develop conditions of optimality on \x). 

The dual problem associated with the SDP of (|l2l ) and (^) (ll| is the problem of maximizing 

D(Z) = -Tv(F Z) (15) 

subject to 

Tv(FiZ) = a, l<i<m; (16) 

Z > 0, (17) 

2 Although typically in the literature the matrices Fi are restricted to be real and symmetric, the SDP formulation 
can be easily extended to include Hermitian matrices Fi\ see e.g., pq . In addition, many of the standard software 
packages for efficiently solving SDPs, for example the Self-Dual-Minimization (SeDuMi) package J2!^, |3(J, allow for 
Hermitian matrices. 
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where Z E B n . A matrix Z E B n is said to be dual feasible if it satisfies (^) and ( |i~7| ) and is strictly 
dual feasible if it satisfies (pi]) and Z > 0. If there exists a strictly feasible point, then the dual 
problem is said to be strictly feasible. 

For any feasible \x) and Z we have that 

m 

P(x) - D(Z) = (c\x) + Tr(F Z) = ^Tr(F^) + Tr(F Z) = Tr(F(x)Z) > 0, (18) 



so that as required, D(Z) < P(x). Furthermore, it can be shown [15] that if both the primal 
problem and the dual problem are strictly feasible, then P = D and \x) is an optimal primal point 
if and only if \x) is primal feasible, and there exists a dual feasible Z E B n such that 

ZF(x) = 0. (19) 

Equation (19) together with (|l6|), ( |T7| ) and ( |l3| ) constitute a set of necessary and sufficient conditions 



for \x) to be an optimal solution to the problem of ( |12|) and (|13|), when both the primal and the 
dual are strictly feasible. 

If Z maximizes D(Z) so that D{Z) = D, then |x) is optimal if and only if F(x) > and 
ZF(x) = 0. 
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3.2 SDP Formulation of Unambiguous Discrimination 



We now show that the unambiguous discrimination problem of ( |To| ) and (11) can be formulated as 
an SDP. Denote by \p) the vector of components pi and by |c) the vector of components — rfc. Then 
our problem is to minimize 

Pip) = (c\p), (20) 



subject to 



1=1 



Pi > 0, 1 < i < m. 



(21) 



To formulate this problem as an SDP, let Fi, < i < m be the block diagonal matrices defined by 



-Qi 



i ■ ■ ■ i L m 



-Qr 



(22) 



Then 



F(p) = F + J2piFi 



i=l 



Xyj=l PiQi 



Pi 



(23) 



so that the constraint F(p) > is equivalent to Yji=iPiQi — ^ r an( ^ Pi > 0, 1 < i < m. Thus the 
problem of ([!(]) and ( |TTj ) reduces to the SDP 



min (c\p) subject to F(p) > 0, 

p£lZ m 



(24) 
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where |c) is the vector of components — r\i with rji being the prior probability of \4>i), and F(p) is 
given by (||). 

To derive a set of necessary and sufficient conditions for optimality on \p), we use the dual 
problem formulation of a general SDP (|l5|)-(|l7|) to formulate the dual problem associated with 



(24), which reduces to 



max -Tr(X), (25) 
xeB r 



subject to 



Ti(QiX)-Zi = 7]i, l<i<m; (26) 
X > 0; (27) 
Zi > 0, 1 < i < m. (28) 

We can immediately verify that both the primal and the dual problem are strictly feasible. 
Therefore it follows that \p) is optimal if and only if the components pi of \p) satisfy (^), there 
exists a matrix X and scalars Zi, 1 < i < m that satisfy (p6|)-(p^), and 



x(i r -f2piQi) = o; (29) 



i=i 



ZiPi = 0, 1 < i < m. (30) 



Note that ( |29| ) implies that for the optimal choice of pi, the largest eigenvalue of YllLiPiQi must 
be equal to 1. This condition has already been derived in 

If X and zi maximize (p5|) subject to (26)-(^8|), then the optimal values of pi can be found by 



solving (p9|) and (^0|) with X = X, Zi = Zi. 
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We summarize our results in the following theorem: 

Theorem 1. Let < i < m} denote a set of state vectors with prior probabilities < 

i < m} in an r- dimensional Hilbert space TC that span an m- dimensional subspace IA of 7i, let 
< i < m } denote the reciprocal states inlA defined by {<pi\(t)k) = b~ik, and let Qi = \4>i) {(f>i\ . 
Let A denote the set of all ordered sets of constants {pi, 1 < i < m} that satisfy pi > and 
YllLiPiQi < Ir> an d let T denote the set of r x r Hermitian matrices X satisfying X > and 
scalars zi > 0, 1 < i < m such that Tr(QiX) — Z{ = rji. Consider the problem min Pig A-P(p) where 
P{p) = — SI^i ViPi an d the dual problem maxjf i2j gr D(X) where D(X) = — Tr{X). Then 

1. For any Pi £ A and X, Zi G T, P(p) > D(X); 

2. There is an optimal \p), denoted \p) , such that P = P{p) < P(p) for any \p) G A; 

3. There is an optimal X and optimal Z{, denoted X and zi, such that D = D(X) > D(X) for 
any X, z; L G T; 

4- P = D; 

5. A set of necessary and sufficient conditions on \p) to minimize P(p) is that pi G A and there 
exists X,Zi G r such that X(I r — Y^iLiPiQi) = and ZiPi = 0, 1 < i < m. 

6. Given X and Z{ a set of necessary and sufficient conditions on \p) to minimize P(p) is that 
Pi G A, X(I r - YliLiPiQi) = an d ZiPi = 0, 1 < i < m. 

As we indicated at the outset, the necessary and sufficient conditions given by Theorem |l] are 
in general hard to solve directly, although they can be used to verify a solution. In addition, these 
conditions can be used to gain insight into the optimal measurement operators. In the next section 
we will use Theorem || to develop necessary and sufficient conditions on a set of state vectors so 
that the EPM is optimal. Contrary to the conditions given by Theorem |], these conditions can be 
easily verified. 
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4 Equal-Probability Measurement 
4.1 Equal-Probability Measurement 

A simple measurement that has been employed for unambiguous state discrimination is the mea- 
surement in which pi = p,l < i < m. This measurement results in equal probability of correctly 
detecting each of the states. We therefore refer to this measurement as the Equal-Probability 
Measurement (EPM). 



To determine the value of p, let <J> have a singular value decomposition (SVD) |22| , 19] of the 
form <I> = UTiV* where U is an r x r unitary matrix, £ is a diagonal r x m matrix with diagonal 
elements <7j > arranged in descending order so that o\ > 02 > • • • > cr m , and V is an m x m 
unitary matrix. Then from @ it follows that 

$ = u(j$yv*, (31) 

where £t is a diagonal m x r matrix with diagonal elements 1/cj. Thus, 

m m 

X>=£ i^x^i = = ^( st )*^*, (32) 
i=i i=i 

and the largest eigenvalue of Y^Li Qi is equal to l/o"^. To satisfy the condition (|29|) the largest 
eigenvalue of pY2i Qi m ust be equal to 1, so that 

P = °m- (33) 



Therefore, our problem reduces to finding necessary and sufficient conditions on the vectors \4>i) 
such that Ilj = cr^Qi minimizes the probability of an inconclusive result. 

In the next section we develop conditions under which the EPM is optimal for unambiguous 
discrimination. In our development, we consider separately the case in which a m has multiplicity 1 
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and the case in which a m has multiplicity greater than 1. We derive a set of necessary and sufficient 
conditions for optimality of the EPM in the first case, and sufficient conditions for optimality in the 
second case. Two broad classes of state sets that satisfy these conditions are discussed in Sections |6] 
and [?]. 

4.2 Conditions For Optimality 

4.2.1 Necessary and Sufficient Conditions 

Let s denote the multiplicity of a m so that a m = cr m _i = • • • = a m+s -i. We first consider the case 
in which s = 1. In this case to satisfy (p9|) and (|27|) we must have that 

X = b\u m )(u m \, (34) 



where \uk) are the columns of U and b > 0. In addition, since pi = p > 0, it follows from ( p0[ ) that 
z i = 0, 1 < i < m so that from pq), 

Tr(QiX) = b\(4>i\u m )\ 2 = m, l<i<m. (35) 

Now, from (|3l|) we have that 

\4>i) = U(^)*\vi), (36) 



where \vi) denotes the ith column of V*. Substituting into (35) 



-jz\vi(m)f =r)i, 1 < i < m, (37) 
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where Vi(k) denotes the A;th component of \vi). Since 



m m 



(38) 



1=1 8=1 



b must be equal to <t 2 



We conclude that when the multiplicity of a m is equal to 1, the EPM is optimal if and only if 
|i>i(m)| 2 = r]i, 1 < i < m, i.e., if and only if each of the elements in the last row of V* is equal to 
the prior probability of the corresponding state. 

4.2.2 Sufficient Conditions 

We now consider the case in which s > 1. To derive a set of sufficient conditions for the EPM to 
be optimal we construct a matrix X that satisfies the conditions of Theorem |]. 



To satisfy (||) and (g7|) we let 




(39) 



k=l 



with bk > 0. Since pi = p > 0, it follows from ( j30|) that z\ = 0, 1 < i < m so that from (|26|), X 



must satisfy 



s 




1 < i < m. 



(40) 



k=l 



Substituting |&) = C7(St)*| Wi ) into @, we have that the constants b^ must satisfy 




1 < i < m 



(41) 



where Vi(k) denotes the A;th component of \vi). 
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We conclude that the EPM is optimal if there exists constants bi > 0,1 < i < s such that 



vi(m) 


I 2 


\vi{m — 


1) 


I 2 ■■ 


\v\{m — 


s + 1) 


I 2 








m 




















h 






\v2{m) 


I 2 


\v2 (m — 


1) 


I 2 ■■ 


\v2(m — 


s + 1) 


I 2 








m 




















b s 






v m (m) 


I 2 


\v m (m - 


i; 


I 2 •• 


■ \v m (m - 


s + r 


l 2 _ 











(42) 



The problem of determining whether there exists a \b) with components bi > such that ( |42| ) 
is satisfied is equivalent to verifying whether a standard linear program is feasible. Specifically, in 
a linear program the objective is to minimize a linear functional of the vector \b) of the form (d\b) 
for some vector \d), subject to the constraints A\b) = \y) andQ \b) > for some given matrix A and 
vector \y). A linear program is feasible if there exists a vector \b) that satisfies the constraints 
Thus we can use standard linear programming techniques to determine whether a \b) exists that 
satisfies (^2|), or equivalently, whether given a set of state vectors with given prior probabilities, 
the EPM is optimal. 

Note, that given a set of state vectors, we can always choose the prior probabilities rji so that 
the EPM is optimal. This follows from the fact that the matrix in (|42"|) depends only on the state 
vectors. Thus, any set of coefficients bi > will give a set of rji > that satisfy (^2|). The coefficients 
r\i will correspond to probabilities if YLi Vi = 1- Since YliLi 1^(^)1 2 = 1 f° r an Vi = Ylt=i ^i, 

and any set of coefficients bi > such that Yli bi = 1 will result in a set of probabilities ryj for which 
the EPM is optimal. 



In [13] the authors raise the question of whether or not cyclic state sets with equal prior 
probabilities are the only state sets for which the EPM is optimal. Here we have shown that the 
EPM can be optimal for any state set, as long as we choose the prior probabilities correctly. In 
Sections || and ^ we consider state sets with equal prior probabilities for which the EPM is optimal, 



The inequality is to be understood as a component-wise inequality. 
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generalizing the result in [[13|], 

We summarize our results regarding the EPM in the following theorem: 

Theorem 2. Let {\(f>i),l < i < m} denote a set of state vectors with prior probabilities {rji,l < 
i < m} in a Hilbert space 7i that span an m- dimensional subspace IA ofH, let {|<^),1 < i < m} 
denote the reciprocal vectors in ti defined by {4>i\<fik) = &ik> an d ^ Q% = I Let $ = UHV* 

denote the matrix of columns \4>i), let \vj) denote the columns ofV* and Vi(k) the kth component 
of \vi), let o~\ > . . . > o m denote the singular values of and let s be the multiplicity of a m . Let 
Hi = cr^Qi denote the equal-probability measurement (EPM) operators. Then, 

1. If s = 1 then the EPM minimizes the probability of an inconclusive result if and only if 
\vi(m)\ 2 = 7]i for 1 < i < m; 

2. If s > 1 then the EPM minimizes the probability of an inconclusive result if there exists 



constants hi > 0, 1 < i < s such that is satisfied; 



3. Given a set of state vectors, we can always choose the prior probabilities rji so that the EPM 
is optimal. Specifically, rji is given by fcj!^) where bi are arbitrary coefficients satisfying b{ > 0, 

and YT=\ b i = !■ 

Theorem || provides necessary and sufficient conditions in the case s = 1 and sufficient conditions 
in the case s > 1 for the EPM to be optimal, which depend on the SVD of $ and the prior 
probabilities rji. It may also be useful to have a criterion which depends explicitly on the given 
states \<fii) and the prior probabilities. Theorem ^ below provides a set of sufficient conditions on 
the states \4>i) and the prior probabilities r]i so that the EPM is optimal. The proof of the Theorem 
is given in the Appendix. In Sections || and [7| we discuss some general classes of state sets that 
satisfy these conditions. 
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Theorem 3. Let {\4>i), 1 < i < m} denote a set of state vectors with prior probabilities {rji, 1 < i < 
m} in a Hilbert space TC that span an m- dimensional subspace U ofTC. Let $ denote the matrix of 
columns \4>i), and let q denote the number of distinct singular values o/<3?. Then the equal-probability 
measurement minimizes the probability of an inconclusive result if (0j|(<l?<J>*)*/ 2-1 |^>j) = r/^at for 
1 < i < m and 1 <t < q, for some constants at- 



5 Computational Aspects 



In the general case there is no closed-form analytical solution to the maximization problem (|20| ) 
subject to (f2l|). However, since this problem is a convex optimization problem, there are very 
efficient methods for solving (f20|). In particular, the optimal vector \p) can be computed on Matlab 
using the linear matrix inequality (LMI) Toolbox. Convenient interfaces for using the LMI toolbox 
are the Matlab packages IQC/? |32| and Self-Dual- Minimization (SeDuMi) p9|, |3(i|] . These algorithms 
are guaranteed to converge to the global optimum in polynomial time within any desired accuracy. 

The number of operations required for each iteration of a general SDP where \x) E lZ m and 
Fi E B n is 0(m 2 n 2 ). However, the computational load can be reduced substantially by exploiting 
structure in the matrices F{. In our problem, these matrices are block diagonal, so that each 
iteration requires on the order of 0(m 4 ) operations [15]. 

To illustrate the computational steps involved in computing the optimal measurement, we now 
consider a specific example. 

Consider the case in which the ensemble consists of 3 state vectors with equal probability 1/3, 
where 



1 

7s 



i 

V2 



1 

V2 





1 
1 



(43) 



To find the optimal measurement operators, we first find the reciprocal states \<pi) . With denoting 
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the matrix of columns \4>i), we have 



1.73 -1.41 

-1.73 1.41 1.41 
1.73 -1.41 



(44) 



and the vectors \cfri) are the columns of <3?. Next, we form the matrices Qi = \4>i){^>i\ which results 



m 





1 


-1 


1 















2 


-2 





Qi = 3 


-1 


1 


-1 


, Q2 = 





2 


-2 


, Qs = 


-2 


2 







1 


-1 


1 







-2 


2 














(45) 



We can now find the optimal vector \p) using the IQC/? package on Matlab. To this end we first 
define the matrices F{ according to (]22]), and define 



(46) 



We then generate the following code, assuming that the matrices Fi and the vector |c) have already 
been defined in Matlab. 
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>> abst_init_lmi 

» p = rectangular (3 , 1) ; 

» F = FO; 

» for i=l:3, 

» evalQ'W = F' num2str(i)]); 
» F = F + p(i)*W; 

» end 

» F > 0; 

>> lmi_mincx_tbx(c ' *p) ; 

» P=value(p) 

The optimal vector \p) is given by 



Vo Initializing the LMI toolbox 
Yo Defining a vector \p) of length 3 
Yo Defining the matrix F(p); here Fi 



Vo Imposing the constraint 

Vo Minimizing (c\p) subject to the constraint 

Vo Getting the optimal value of p 




0.17 
0.17 



(47) 



and the optimal measurement operators IT = piQi are 



























0.34 


-0.34 





n x = 











, n 2 = 





0.34 


-0.34 


, n 3 = 


-0.34 


0.34 





















-0.34 


0.34 














(48) 



We can now use the necessary and sufficient conditions derived in Section 3.2 and summarized 
in Theorem |l] to verify that \p) given by (|47| ) is the optimal probability vector. To this end we first 
form the matrix T = I r — X/f=i n»- Using the eigendecomposition of T we conclude that the null 
space of T has dimension 1 and is spanned by the vector 



u 



-0.81 
0.41 
0.41 



(49) 



Therefore to satisfy ( |29| ) and (|27|), X must be equal to X = a\u)(u\ for some a > 0. Since p\ = 



and P2,P3 > 0, (30) and (|28| ) imply that Z2 = £3 = and Z\ > 0. Therefore, from (|2q) we must 
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have that 



and 



Tr(Q 2 X) = Tr(Q 3 X) 



3' 



To satisfy ( |50|) we choose 



1 



0.11. 



(50) 



(51) 



(52) 



3(u\Q 2 \u) 

With this choice of a, Tr(Q3X) = 1/3 and Tr(QiX) = 0.89 > 1/3, so that the necessary and 
sufficient conditions are satisfied. 

Now, suppose that instead of equal prior probabilities we assume that the prior probabilities 
are r\\ = 0.6, 772 = 0.2, 773 = 0.2. These priors where chosen to be equal the elements of the last row 
of V* where $ = UTiV*. Since the smallest square singular value of <£, cr| = 0.07, has multiplicity 



1, (42) is satisfied and the EPM consisting of the measurement operators IT = pQi with p = 0.07, 
minimizes the probability of an inconclusive result. As before, we can immediately verify that this 
is indeed the correct solution using the necessary and sufficient conditions of Theorem |l]. For this 
choice of IT, T = I r — p ^4=1 Qi> an d the null space of T is spanned by the vector 



u 



0.68 
-0.52 
-0.52 



(53) 



Therefore X must be equal to X = a\u)(u\ for some a > 0. Since pi = p > for all i, Z{ = 0, 1 < 
i < 3 so that we must have 



TV(QiX) = 0.6, Tr(Q 2 X) = 0.2, Tr(Q 3 X) = 0.2. 



(54) 
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If we choose a = 0.6/ (u\Qi\u) = 0.07, then ( p4[ ) is satisfied, and the EPM is optimal. 

In the remainder of the paper we use the sufficient conditions of Theorem || to derive the 
optimal unambiguous measurement for state sets with certain symmetry properties. The symmetry 
properties we consider are quite general, and include many cases of practical interest. Specifically, 
in Section ^ we consider geometrically uniform state sets, and in Section we consider compound 
geometrically uniform state sets. It is interesting to note that for these classes of state sets, the 
optimal measurement that minimizes the probability of a detection error is also known explicitly 

MM- 

6 Geometrically Uniform State Sets 

In this section we consider the case in which the state vectors \(f>i) are defined over a group of unitary 
matrices and are generated by a single generating vector. Such a state set is called geometrically 
uniform ( GU) [plf . We first obtain a convenient characterization of the EPM in this case and then 
show that the EPM is optimal. This result generalizes a similar result of Chefles and Barnett [O], 

6.1 GU State Sets 

Let Q be a finite group of m unitary matrices Ui on 7i. That is, Q contains the identity matrix 
I r ; if Q contains Ui, then it also contains its inverse Uf = U*; and the product UUj of any two 
elements of Q is in Q [33|. 

A state set generated by Q using a single generating vector \<j>) is a set S = {]&) = Ui\cj)),Ui G G}- 
The group Q will be called the generating group of S. For concreteness we assume that U\ = I r so 
that \4>i) = \4>). Such a state set has strong symmetry properties and is called GU. For consistency 
with the symmetry of S, we will assume equiprobable prior probabilities on S. 

Alternatively, a state set is GU if given any two states \4>i) and \<j>j) in the set, there is an 
isometry (a norm-preserving linear transformation) that transforms \4>i) into while leaving the 
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set invariant [18|. Intuitively, a state set is GU if it "looks the same" geometrically from any of the 



states in the set. Some examples of GU state sets are considered in [18, IS]. 

We note that in |l9| a GU state set was defined over an abelian group of unitary matrices. Here 
we are not requiring the group Q to be abelian. 

A cyclic state set is a special case of a GU state set in which the generating group Q has elements 
Ui = 1 < i < m, where Z is a unitary matrix with Z m = I r . A cyclic group generates a cyclic 
state set S = {\4>i) = Z t ~ 1 \(j)), 1 < i < m}, where \4>) is arbitrary. 

Any binary state set S = {|</>i), \ <p2)} is a GU cyclic state set, because it can be generated by 
the binary group Q = {7 r ,-R}, where R is the reflection about the hyperplane halfway between the 
two states. Since R represents a reflection, R is unitary and R 2 = I r . 

6.2 The EPM for GU States 

To derive the EPM for a GU state set with generating group Q, we need to determine the reciprocal 



states \4>i). It was shown in [21, |2(| that for a GU state set with generating group G, commutes 



with each of the matrices Ui £ Q. For completeness we repeat the argument here. Expressing 



as 



8=1 i=l 



we have that for all j, 



i=l 

m 



8=1 



8=1 



Uj^*, (56) 
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since {U*Ui, 1 < i < m} is just a permutation of Q. 

If commutes with Uj, then T = ($$*)^ also commutes with Uj for all j. Thus, from (||) 
the reciprocal states are 

|&) = = TUi\<f>) = UiT\(t>) = Ui\4>), (57) 

where 

|0)=T|0) = ($$*)t|0). (58) 

It follows that the reciprocal states are also GU with generating group Q and generating vector 
given by (p8|). Therefore, to compute the reciprocal states for a GU state set all we need is to 
compute the generating vector | <ft) . The remaining vectors are then obtained by applying the group 
Q to \4>). The EPM is then given by the measurement operators 

Q i =pU i \$){0\U i , (59) 

where p is equal to the smallest eigenvalue of $$*. 
6.3 Optimality of the EPM 

We now show that the EPM is optimal for GU state sets with equal prior probabilities r\i = 1/m. 
Since commutes with Uj for all j, (<]?<3?*) a also commutes with Uj for any a. Therefore for all 
t, 

Since (^iK^^*)*/ 2-1 !^) does not depend on i, from Theorem || we conclude that the EPM is 
optimal. 
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We summarize our results regarding GU state sets in the following theorem: 

Theorem 4 (GU state sets). Let S = {\4>i} = Ui\4>),Ui € Q} be a geometrically uniform (GU) 
state set generated by a finite group Q of unitary matrices, where \<p) is an arbitrary state, and 
let <& be the matrix of columns \4>i) . Then the measurement that minimizes the probability of an 
inconclusive result is equal to the equal-probability measurement, and consists of the measurement 
operators 

H =p\4>i)(4>i\, 



where {\4>i) = Ui\<fi), Ui £ Q}, 



\4>) = (***?\4>), 



and p is the smallest eigenvalue o/$$*. 

6.4 Example Of A GU State Set 

We now consider an example of a GU state set. 

Consider the group Q of m = 4 unitary matrices Ui, where 



U!=h, U 2 



1 









-1 

1 









-1 



, U 3 



1 

1 










-1 









-1 



, u 4 = u 2 u 3 . 



(61) 
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Let the state set be S = {|<&) = Ui\(f>), 1 < % < 4}, where |0) = l/(3\/2)[2 2 1 3]*, so that 



3\/2 



2 
2 
1 
3 



2 2 

-2 2 

1 -1 

-3 -3 



(62) 



From Theorem |] the measurement that minimizes the probability of an inconclusive result is the 
EPM. Furthermore, the reciprocal states \4>i) are also GU with generating group Q and generator 



($$*)t| 



4^/2 



3 
3 
6 
2 



(63) 



so that {\4>i) = Ui\4>), 1 < i < 4}. Since 



<j><j)* = _ 
9 



4 

4 

10 

9 



(64) 



p = 2/9 and the EPM measurement operators are Ik = (2/9)Qi = (2/9)Ui\4>){4>\U*. 

We can now use the necessary and sufficient conditions of Theorem [l] to verify that IT = 
(2/9)\<f>i){4>i\ are indeed the optimal measurement operators. To this end we first form the matrix 
T = I r — ^2f = i Ilj. Using the eigendecomposition of T we conclude that the null space of T has 
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dimension 1 and is spanned by the vector 



(65) 



Therefore to satisfy (^) and (p7|), X must be equal to X = a\u){u\ for some a > 0. Since 
Pi = 2/9 > 0, 1 < i < 4, @) and (||) imply that zi = 0, 1 < i < 4. Therefore, from @ we must 
have that 



Tr(QiX) = Tr(Q 2 X) = Tr(Q 3 X) = Tr(Q 4 ^) 



1 



To satisfy (|66|) we choose 



1 



(66) 



(67) 



A{u\Q x \u) 9 

With this choice of a, Tr(Q2-^") = Tr(Q3^) = ^(Q^X) = 1/4, so that as we expect the necessary 



and sufficient conditions are satisfied. 



7 Compound Geometrically Uniform State Sets 

In Section [6| we showed that the optimal measurement for a GU state set is the EPM associated 
with this set. We also showed that the reciprocal states are themselves GU and can therefore be 
computed using a single generator. In this section, we consider state sets which consist of subsets 
that are GU, and are therefore referred to as compound geometrically uniform (CGU) [^TJ. As we 
show, the reciprocal states are also CGU so that they can be computed using a set of generators. 
Under a certain condition on the generating vectors, we also show that the EPM associated with a 
CGU state set is optimal. 

A CGU state set is defined as a set of vectors S = {\4>ik),l < i < 1,1 < k < r} such that 
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\<t>ik) = Ui\(f>k), where the matrices {Ui, 1 < i < 1} are unitary and form a group G, and the vectors 
{|0ifc}) 1 < k < r} are the generating vectors. For consistency with the symmetry of S, we will 
assume equiprobable prior probabilities on S. 

A CGU state set is in general not GU. However, for every k, the vectors {\4>ik), 1 < i < 1} are 
a GU state set with generating group Q. Examples of CGU state sets are considered in [pi], |2Q| . 

7.1 The EPM for CGU State Sets 

We now derive the EPM for a CGU state set with equal prior probabilities. Let <3? denote the 
matrix of columns \4>ik), where the first I columns correspond to k = 1, and so forth. Then for a 
CGU state set with generating group Q, it was shown in [21, that commutes with each of 



the matrices Ui G Q. If commutes with Ui, then T = ($$*)^ also commutes with Ui for all i. 
Thus, the reciprocal states are 

\4>ik) = T\<j> ik ) = TUi\<f> k ) = UiT\<f> k ) = Ui\$ k ), (68) 



where 

= T|0 fc ) = (69) 

Therefore the reciprocal states are also CGU with generating group Q and generating vectors \4>k) 
given by (|69|). To compute these vectors all we need is to compute the generating vectors \4>k)- The 
remaining vectors are then obtained by applying the group Q to each of the generating vectors. 

7.2 CGU State Sets With GU Generators 

A special class of CGU state sets is CGU state sets with GU generators plj in which the generating 
vectors {|0fc),l < k < r} are themselves GU. Specifically, {\4>k} = ^fcl^)} f° r some generator \<p), 
where the matrices {V^, 1 < k < r} are unitary, and form a group Q. Examples of CGU state sets 
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with GU generators are considered in [20]. 

Suppose that U{ and V k commute up to a phase factor for all i and k so that UiV k = V k Uie^^' k ^ 
where 9(i, k) is an arbitrary phase function that may depend on the indices i and k. In this case 
we say that Q and Q commute up to a phase factor (in the special case in which 9 = so that 
UiVk = VkUi for all i, k, the resulting state set is GU plf ). Then for all i,k, <I )C I ) * commutes with 
UiVk 20]. The reciprocal states \<j>ik) of the vectors \4>ik) are therefore given by 



TUiV k \4>) = UiV k T\4>) = UiV k \4>), (70) 



where |0) = T\(f>). Thus even though the state set is not in general GU, the reciprocal states can 
be computed using a single generating vector. 

Alternatively, we can express \4>ik) as |<^ifc) = Ui\(f> k ) where the generators \(j) k ) are given by 

\h) = v k \4>). (7i) 

From (71) it follows that the generators \(p k ) are GU with generating group Q = {V k , 1 < k < r} 



and generator 

We conclude that for a CGU state set with commuting GU generators and generating group Q, 
the reciprocal states are also CGU with commuting GU generators and generating group Q. 

7.3 The Optimal Measurement for CGU State Sets Satisfying a Weighted Norm 
Constraint 

We now show that if the generating vectors \cf>k) satisfy 

(McD^)*/ 2 - 1 !^) =<H, l<k<r,l<t<q, (72) 
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where q is the number of distinct singular values of <3?, then the EPM is optimal. 
Prom Theorem |3] it follows that it is sufficient to show that (ff^ ) implies 



|(***)*/ 2 - 1 |0 tt )=a t , l<i<l,l<k<r,l<t<q. (73) 



Now, 

($*$)*/ 2 - 1 |0 ifc ) = (M)*/ 2 " 1 ^) = (74) 

so that 

(^K*^)*/ 2 - 1 !^) = (^it^^C***)*/ 2 - 1 !^) = (^|($*$) t / 2 - 1 |0fc) = <h, (75) 



establishing (|73|). 

For CGU state sets with GU generators {\4>k) = V^. 1 0} } where Vk £ Q and Q and Q commute 
up to a phase factor, the EPM is optimal. This follows from the fact that in this case (^) is always 



satisfied. To see this, we first note that commutes with for each k [21]. Therefore for all k, 



* fc |(***)*/2-%) = (0|y fe *($**)*/ 2 -V fe |^) = (^I^VfcC***)*/ 2 - 1 ^) = (0|($*$) t/2_1 |0>. (76) 
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We summarize our results regarding CGU state sets in the following theorem: 

Theorem 5 (CGU state sets). Let S = {\<j>ik) = Ui\(j>k), 1 < % < 1, 1 < k < r}, be a compound 
geometrically uniform (CGU) state set generated by a finite group Q = {Ui, 1 < i < 1} of unitary 
matrices and generating vectors {\4>k), 1 < k < r}, and let $ be the matrix of columns \4>ik)- Then 
the equal-probability measurement (EPM) consists of the measurement operators 

n» = p\4>ik){4>ik\, 

where {\4>ik) = Ui\<t> k ), 1 < i < 1, 1 < k < r}, 

and p is equal to the smallest eigenvalue of . 
The EPM has the following properties: 

1. If {<fik\(&&*) t ^ 2 ~ 1 \4 > k) = a t for 1 < k < r, 1 < t < q where q is the number of distinct 
eigenvalues of , then the EPM minimizes the probability of an inconclusive result. 

2. If the generating vectors {\<ftk) = VJt | ^6) , 1 < k < r} are geometrically uniform with UiVk = 
V k Uie? e ^ for alli,k, then 

(a) \4>ik) = UiVk\4>) where \cj>) = (<E><I>* y\4>) so that the reciprocal states are CGU with geo- 
metrically uniform generators; 

(b) The EPM is optimal; 

(c) If in addition 6(i, k) = for all i, k, then the vectors {(pik, 1 < i < Z,l < k < r} form a 
geometrically uniform state set. 
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Appendix 

Proof of Theorem |3] 

In this appendix we prove Theorem |j| 

Let Aj, 1 < i < q denote the singular values of $ without multiplicity so that Ai = o\ and 
A 9 = cr m , and let Si denote the multiplicity of Aj. Define 



and 



H 



E£iM*)|- 



A 



Ai A2 
\\ A| 



A' A2 



A„ 



\2 

A <2 



Ei=i Pi\vi(m - s q + i)\ 2 Y,iLx Pi\v2{m - s q + 



Ei=i M*i +■ 



Eiii A Mm - Sg + i) 



(77) 



(78) 



for some > 0. Finally, let N be the matrix with ith column equal to r]i\a) where \a) is an 
arbitrary vector. 

Now, suppose that AH = N. Then A\hi) = rji\a) where |/tj) denotes the ith. column of H. Since 
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A is invertible, this implies that 



— hi(k) = — hj(k) 1 < i. j < m, 1 < k < q. 



(79) 



For k = q, ( |79D reduces to We therefore conclude that a sufficient condition for the EPM to 

be optimal is that AH = N for some > 0. Taking = 1 for each i, we can express ^4ff as 



AH 



<J\ &2 ■■■ Cm 

2 



CTr, 



Ml)| 2 \v 2 (l)\ 2 ■■■ \v 2 (m)\ 2 
h(2)| 2 M2)| 2 ••• \v 2 (2)\ 2 

\vi (m)\ 2 |f2("i-)| 2 ••• b2(w)| 2 



Y. 



Then we have that 



Y U = j2^Mi)\ 2 = m^*) t/2 - 1 \^i). 



1=1 



Therefore AH = N reduces to the condition that 



(80) 



(81) 



^Kfcfc*)*/ 2 - 1 !^) = mat, l<l<m,l<t<q, 



(82) 



for some constants at- 



34 



References 

[1] C. H. Bennett and P. W. Shor, "Quantum information theory," IEEE Trans. Inform. Theory, 
vol. 44, no. 6, pp. 2724-2742, 1998. 

[2] A. Peres, Quantum Theory: Concepts and Methods, Boston: Kluwer, 1995. 

[3] C. W. Helstrom, Quantum Detection and Estimation Theory, New York: Academic Press, 
1976. 

[4] A. S. Holevo, "Statistical decisions in quantum theory," J. Multivar. Anal, vol. 3, pp. 337-394, 
Dec. 1973. 

[5] H. P. Yuen, R. S. Kennedy, and M. Lax, "Optimum testing of multiple hypotheses in quantum 
detection theory," IEEE Trans. Inform. Theory, vol. IT-21, pp. 125-134, Mar. 1975. 

[6] Y. C. Eldar, A. Megretski, and G. C. Verghese, "Designing optimal quantum detectors 
via semidefinite programming," IEEE Trans. Inform. Theory, to appear; also available at 
http : //www . arXiv . org/abs/ quant -ph/ 0205178. 

[7] I. D. Ivanovic, "How to differentiate between non-orthogonal states," Phys. Lett. A, vol. 123, 
pp. 257-259, Aug. 1987. 

[8] D. Dieks, "Overlap and distinguishability of quantum states," Phys. Lett. A, vol. 126, pp. 
303-307, 1988. 

[9] A. Peres, "How to differentiate between non-orthogonal states," Phys. Lett. A, vol. 128, pp. 
19, Mar. 1988. 

[10] G. Jaeger and A. Shimony, "Optimal distinction between two non-orthogonal quantum states," 
Phys. Lett. A, vol. 197, pp. 83-87, 1995. 

[11] A. Peres and D. R. Terno, "Optimal distinction between non-orthogonal quantum states," J. 
Phys. A, vol. 31, pp. 7105-7111, 1998. 

[12] A. Chefles, "Unambiguous discrimination between linearly independent quantum states," 
Phys. Lett. A, vol. 239, pp. 339-347, Apr. 1998. 

[13] A. Chefles and S. M. Barnett, "Optimum unambiguous discrimination between linearly inde- 
pendent symmetric states," Phys. Lett. A, vol. 250, pp. 223-229, 1998. 

[14] A. Peres, "Neumark's theorem and quantum inseparability," Found. Phys., vol. 20, no. 12, 
pp. 1441-1453, 1990. 

[15] L. Vandenberghe and S. Boyd, "Semidefinite programming," SIAM Rev., vol. 38, no. 1, pp. 
40-95, Mar. 1996. 

[16] Y. Nesterov and A. Nemirovski, Interior- Point Polynomial Algorithms in Convex Program- 
ming, Philadelphia, PE: SIAM, 1994. 



35 



[17] F. Alizadeh, Combinatorial Optimization With Interior Point Methods and Semi-Definite 
Matrices, Ph.D. thesis, University of Minnesota, Minneapolis, MN, Oct. 1991. 

[18] G. D. Forney, Jr., "Geometrically uniform codes," IEEE Trans. Inform. Theory, vol. 37, pp. 
1241-1260, Sep. 1991. 

[19] Y. C. Eldar and G. D. Forney, Jr., "On quantum detection and the square-root measurement," 
IEEE Trans. Inform. Theory, vol. 47, pp. 858-872, Mar. 2001. 

[20] Y. C. Eldar, A. Megretski, and G. C. Verghese, "Optimal detection of symmetric mixed 
quantum states," in preparation. 

[21] Y. C. Eldar and H. Bolcskei, "Geometrically uniform frames," IEEE Trans. Inform. Theory, 
to appear; also available at http://arXiv.org/abs/math.FA/0108096. 

[22] G. H. Golub and C. F. Van Loan, Matrix Computations, Baltimore MD: Johns Hopkins Univ. 
Press, third edition, 1996. 

[23] M. Jezek, J. Rehacek, and J. Fiurasek, "Finding optimal strategies for minimum-error 
quantum-state discrimination," http: //www. arXiv. org/abs/quant-ph/0201109. 

[24] A. C. Doherty, P. A. Parrilo, and F. M. Spedalieri, "Distinguishing separable and entangled 
states," http : //www . arXiv . org/abs/quant-ph/0112007. 

[25] E. M. Rains, "A semidefinite program for distillable entanglement," IEEE Trans. Inform. 
Theory, vol. 47, no. 7, pp. 2921-2933, Nov. 2001. 

[26] K. Audenaert and B. De Moor, "Optimizing completely positive maps using semidefinite 
programming," http : //www . arXiv . org/abs/quant-ph/0109155. 

[27] X. Sun, S. Zhang, Y. Feng, and M. Ying, "Mathematical nature of and a family of lower 
bounds for the success probability of unambiguous discrimination," Phys. Rev. A, vol. 65, 
Sep. 2002. 

[28] M. X. Goemans and D. P. Williamson, "Approximation algorithms for MAX-3-CUT and 
other problems via complex semidefinite programming," in ACM Symposium on Theory of 
Computing, 2001, pp. 443-452. 

[29] J. F. Strum, "Using SeDuMi 1.02, a MATLAB toolbox for optimization over symmetric cones," 
Optimization Methods and Software, vol. 11-12, pp. 625-653, 1999. 

[30] D. Peaucelle, D. Henrion, and Y. Labit, "Users guide for SeDuMi interface 1.03," available at 
http : //www . laas . f r/peaucell/SeDuMilnt .html. 

[31] D. Bertsimas and J. Tsitsiklis, Introduction to Linear Optimization, Belmont MA: Athena 
Scientific, 1997. 

[32] A. Megretski, C-Y. Kao, U. Jonsson, and A. Rantzer, "A guide to IQC/3: Software for 
robustness analysis," available at http://web.mit.edu/cykao/www. 

[33] M. A. Armstrong, Groups and Symmetry, New York: Springer- Verlag, 1988. 



36 



